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Abstract 

In this paper we show the capabihties of the Large Volume Detector (INFN Gran 
Sasso National Laboratory) to identify a neutrino burst associated to a supernova 
explosion, in the absence of an "external trigger", e.g., an optical observation. We 
describe how the detector trigger and event selection have been optimized for this 
purpose, and we detail the algorithm used for the on-line burst recognition. The 
on-line sensitivity of the detector is defined and discussed in terms of supernova 
distance and intensity at the source. 
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1 Introduction 



The detection of neutrinos from SN1987A marked the beginning of a new 
phase of neutrino astrophysics [T][2][5|l] . In spite of the lack of a firmly es- 
tablished theory of core collapse supernova explosion, the correlated neutrino 
emission appears to be well established. However, since this first v observation 
was guided by the optical one, the detector capabilities of identifying a z/ burst 
in the absence of an " external trigger" should be demonstrated very carefully. 
In the presence of an electromagnetic counterpart, on the other hand, the 
prompt identification of the neutrino signal could alert the worldwide network 
of observatories allowing study of all aspects of the rare event from its onset. 
The Large Volume Detector (LVD), in the INFN Gran Sasso National Labo- 
ratory (Italy), at the depth of 3600 m w.e., is a 1 kt liquid scintillator detector 
whose major purpose is monitoring the Galaxy to study neutrino bursts from 
gravitational stellar collapses [S]. Besides interactions with protons and car- 
bon nuclei in the liquid scintillator, LVD is also sensitive to interactions with 
the iron nuclei of the support structure whose total mass is 0.9 kt [6J. The 
experiment has been taking data, under different configurations, since 1992, 
reaching in 2001 its present and final configuration. Its modularity and rock 
overburden, together with the trigger strategy, make this detector particularly 
suited to disentangle on-line a cluster of neutrino signals from the background. 
We will discuss in this paper the LVD performances from the point of view of 
the on-line identification of a neutrino burst: we will describe in section 2 the 
trigger of the detector, the event selection and the on-line burst recognition. In 
section 3 we will define and discuss the detector sensitivity to neutrino bursts 
which, as we will show in section 4, can be expressed in terms of supernova 
distance or neutrino intensity at the source. 



2 The LVD event selection chain 

2.1 The trigger 

LVD consists of an array of 840 scintillator counters, 1.5 m^ each. The whole ar- 
ray is divided in three identical "towers" with independent high voltage power 
supply, trigger and data acquisition (see figured]). In turn, each tower consists 
of 35 "modules" hosting a cluster of 8 counters. Each counter is viewed from 
the top by three 15 cm photomultiplier tubes (PMTs) FEU49b or FEU125. 
The charge of the summed PMT signals is digitized by a 12 bit ADC (conver- 
sion time = 1 /is) and the arrival time is measured with a relative accuracy of 
12.5 ns and an absolute one of 100 ns [7]. 

The main neutrino reaction in LVD is VeV ~^ ^"'"n, which gives two detectable 
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signals: the prompt one due to the e"*" (visible energy E^is ~ Ep^ — 1.8 MeV 
+ 2 nieC^) followed by the signal from the np d'y capture {E^ = 2.2 MeV, 
mean capture time ~ 185 /xs,). The trigger logic is optimized for the detec- 
tion of both products of the inverse beta decay and is based on the three-fold 
coincidence of the PMTs of a single counter. Each PMT is discriminated at 
two different thresholds resulting in two possible levels of coincidence between 
a counter's PMTs: H and L, corresponding to Sh — 4 MeV and Sl — 1 
MeV. The H coincidence in any counter represents the trigger condition for 
the array. Once a trigger has been identified, the charge and time of the three 
summed PMTs' signals are stored in a memory buffer. All signals satisfying 
the L coincidences in the same module of the trigger counter are also stored, if 
they occur within 1 ms. The average efficiency of the H trigger (for electrons) 
is shown in fig. [2] as a function of the visible energy -Evis- The average neutron 
detection efficiency, en, amounts to about 50% for neutrons detected in the 
same counter where the positron has been detected [S]. 




Fig. 1. LVD schematic view: the 
building block is a module of 8 scin- 
tillator counters, 35 modules form 
a tower. 
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Fig. 2. Trigger efficiency, averaged 
over all the scintillator counters, 
for the H coincidence versus visi- 
ble energy ^^is- 



One millisecond after the occurrence of a trigger, the memory buffers ^ I. con- 
taining the charge and time information of both H and L signals, are read out. 
This is performed independently on the three towers, without introducing any 
dead time. The tower event fragments, if present, are then sent to a central 
processor which provides event-building requiring time coincidence of the first 



Each memory buffer, one per module, is able to store up to 512 signals, corre- 
sponding to 50000 signals in the whole apparatus. 
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pulse registered in each tower within 10 microsecond. The whole duration of 
one event can extend up to 1 ms. 

2.2 Event selection criteria 

The basis of the search for neutrino bursts is the identification of clusters 
of signals in fixed time windows. A preliminary step consists in the selection 
of good signals detected by good counters; we thus apply to the events two 
different series of cuts: (A) to the signals, (B) to the counters, namely: 

(A) Cut on signals 

(1) The energy of the signals must be in the range Ecut < Esignai <100 MeV. 
Two energy thresholds are considered: Ecut = 7 MeV, corresponds to 
an average trigger efficiency > 90%, and Ecut = 10 MeV, to an average 
trigger efficiency > 95% (see fig. [21) ■ 

(2) Events with signals in coincidence in 2 or more counters of the array, 
within 200 ns, are rejected since they are considered muon candidates ^ I. 
The probability for electron or positron with energy up to 60 MeV to 
trigger 2 LVD counters is less than 10% [9]. 

(B) Cut on counters 

(1) The response to muons of atmospheric origin is used to identify and 
discard not properly working counters. The measurement of muons as- 
sociated with the CNGS beam [10] provides complementary tests of the 
detection capabilities of the detector. A wrong muon counting rate or 
spectrum can be due to scintillator, PMT or electronics problems (ADC 
or TDC). Counters rejected for this reason are less than 5% of the to- 
tal and represents a steady loss of active mass requiring a maintenance 
intervention. 

(2) Counters with a background rate (for E > 7 MeV) R > 3 • 10"^ s~^, 
during the last two hours of operation, are rejected as noisy. The problem 
arises from electronics or bad energy calibration and usually regards less 
than 2% of the counters. Also in this case counters involved are almost 
always the same, a new energy calibration or a maintenance intervention 
are required. The typical counting rate distribution for all good counters 
(after cuts), averaged over 150 days of data taking, is shown in figure O 

(3) Counters that take part too often in a cluster of signals are rejected, and 
the cluster redefined. Namely, if m is the cluster multiplicity and the 
number of active counters, a counter will be rejected if its multiplicity mi 

^ This rejection corresponds to a dead time ~ 0.01 %, taking into account that the 
global muon rate in LVD is about 0.1 muon s^^. 
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corresponds to a probability I^k^m- P(k;ni/Nc)< 1 • 10~^. This is due to 
sporadical (of the order of once per month) and local electric noise and 
determines a mass loss ^ 0.5%. 

The effect of the cuts on counters is to adjust dynamically the LVD active 
mass, Mact, by rejecting counters that are steadily out of order (1 and 2) 
and/or that are sporadically noisy (3). A snapshot of the detector active mass 
in the last two years is shown in figure IH In the same period the average 
background counting rate of the whole array was fbk=0.2 Hz if Ecut=7 MeV 
and fbk=0.03 Hz for Ecut=10 MeV. 




Trigger Rate (10 Hz) 

Fig. 3. Counting rate distribution 
(for Evis > 7 MeV) for all good 
counters (after cuts), averaged 
over 150 days of data taking. The 
average rate is ^ = 2 • 10^^ s^^ 
counter"^. 
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Fig. 4. On-line LVD active mass 
during the last two year of oper- 
ation. 



At the end of the whole selection procedure the time distribution of the signals 
is well described by the Poisson statistics, as can be seen in figure [5l where we 
show the difference of the arrival time between successive signals, and in figure 
[6] where we show the fluctuations of the 5 min counting rate, S5, with respect to 
the local average value (measured during a time window of 40 min) in units of 
the expected error, a^xp, calculated assuming pure Poisson fluctuations. The 
typical experimental distribution, obtained during 100 days of operation, is 
fitted with a Gaussian one with mean equal zero and a = 1.01, showing that 
the residual non-Poissonian contribution to the fluctuations (Tres = a/o^ — 1 is 



^ For example: in a cluster (m=100, At=20 s) detected when there are 800 active 
counters, a counter will be rejected (and the cluster redefined) if 5 or more signals 
inside the cluster come from it. 
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less than 15%. 
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Fig. 5. Distribution of the difference 
of the arrival time between succes- 
sive signals, compared with Pois- 
son expectations. 
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Mean = 0.00±0.02 
Sigma = 1.01 ±0.01 




Fig. 6. Distribution of the fluctua- 
tions of the 5 minutes counting rate 
( fs = (s5 - S5)/o-oxp)- The superim- 
posed curve is the free parameters 
Gaussian fit. The statistics repre- 
sents 100 days of data taking. 



2.3 The u hursts selection algorithm 



The core of the algorithm for the on-line selection of candidate neutrino bursts 
is the search for a cluster of H signals within a fixed-duration time window, 
At. The candidate burst is simply characterized by its multiplicity m, i.e., the 
number of pulses detected in At, and by At itself. All the other characteristics 
of the cluster, e.g., detailed time structure, energy spectra, z/ flavor content and 
topological distribution of signals inside the detector are left to a subsequent 
independent analysis. Based on this principle, the LVD data are continuously 
analyzed by an on-line "supernova monitor". In detail, each data period, T, 
is scanned through a "sliding window" with duration At = 20 s, that is, it 
is divided into N = 2 --^ — 1 intervals, each one starting in the middle of 
the previous one, so that the unbiased time window is 10 s. The frequency of 
clusters of duration 20 s and multiplicity > m, due to background, is: 

oo /■ 

Fi^{m, /bk, 20 s) = 8640 ■ ^ P{k- 20 ■ ^) event ■ day-' (1) 

k>m 
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where fbk is the background counting rate of the detector for Evis > Ecut, 
Pik; fhk'^t) is the Poisson probabihty to have clusters of multiphcity k if 
/bkAt is the average background multiphcity, and 8640 is the number of trials 
per dayP] For example, a cluster with m=10 can be produced by background 
fluctuations once every 100 years if fbk = 0.03 s^^; to have the same signifi- 
cance with a higher background, e.g. fbk = 0.2 s^^, a cluster with much higher 
multiplicity, m=22, would be required (see figure [7| where m versus is 
shown, for two different background conditions). For these background rates, 
m = 10 and m = 22 correspond to the minimum multiplicity, m^in, to have 
an imitation frequency < 1 ■ 10~^ yr~^. 

In LVD the search for burst candidates is performed for both energy cuts: 
7 and 10 MeV. The chosen Fim, below which the detected cluster will be an 
on-line candidate supernova event, is 1 per 100 year working stand-alone while 
it is relaxed to 1 per month working in coincidence with other detectors, as in 
the SNEWS project [II]. 

In figure [8] we show the distributions of the observed time intervals between 
selected clusters at the imitation frequency = 1/day and = 1/month 
during 688 days (between July 5*^^, 2005 and May 23'^'^, 2007) and E^ut = 7 
MeV. The mean rates, derived from Poisson fits, are F°^^ = 1.24 day^^ (with 
"^min varying between 13 and 15) and = 1.28 month^^ (^min between 
15 and 18), respectively, meaning that, within 25%, the time behavior of the 
background is consistent with expectations. 

We can conclude that the background trend is predictable even during long pe- 
riods of data acquisition and with variable detector conditions, allowing us to 
define the significance of a neutrino burst in terms of imitation frequency, Fim. 



^ Once a candidate cluster (m, 20s) has been identified, the algorithm will search 
for the most probable starting point of the signal within the time window. For all 
possible sub-clusters with multiplicity 2 < k < m and duration < Ai < 20s the 
corresponding Poisson probability Pk>m is calculated and the absolute minimum 
identified. The time of the first event of the least probable sub-cluster is assumed 
as the start time of the signal. 
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Fig. 7. Minimum cluster multiplic- 
ity rrijam vs. the imitation fre- 
quency Fim- Triangles correspond 
to Ecut=7 MeV and fbk=0.2 Hz, 
while squares to Ecut=10 MeV 
and fbk=0.03 Hz. The two verti- 
cal lines represent the thresh- 
olds of 1 candidate per 100 year 
(stand-alone) and of 1 per month 
(SNEWS). 
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Fig. 8. Distribution of the time in- 
tervals between observed clusters 
(histograms) fitted by Poisson laws 
(lines) for = 1/day (solid) and 
^im = 1/month (dashed), during 
688 days and Ecut = 7 MeV. 



3 Sensitivity to supernova neutrino bursts 

In this section we discuss the sensitivity of the described on-hne selection al- 
gorithm to the recognition of a supernova event. The selection method defines 
a candidate as any cluster of m > mmin signals within a window of At = 20 s. 
For a known background rate, mmin corresponds to a chosen which is set 
as a threshold. This multiplicity represents the minimum number of neutrino 
interactions required to produce a supernova " alarm" , and contains two terms, 
one due to the background, /bkAt, and the other due to the neutrino signal. 
In particular, for LVD, considering only inverse beta decay (IBD) reactions, 
which are the dominant ones at least in the "standard" supernova model, and 
simply approximating the detector response to -Evis = Ep^ — 0.8 MeV, one can 
write: 

lOs lOOMeV 

mmin = /bkAt + M,,t iVp e(^cut) J dt J HE^^^t) ■ a{EpJdEp^{2) 

Ecut+0.8MeV 

where: Mact is the active mass, A^p = 9.34 10^^ is the number of free protons 
in a scintillator ton, e(£'cut) is the trigger efficiency approximated as constant 
(e = 0.9 for E^^t = 7 MeV and e = 0.95 for E^^t = 10 MeV, see fig. ED, 
cr{EpJ is the IBD cross section[12] and ^{Ep^,t) the differential z/g intensity 
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at the detector. The upper hmit in the time integral (10 s) corresponds to the 
maximum unbiased cluster duration. 

Hence the integral on the right side of (2) is the detector burst sensitivity, S, 
in terms of minimum neutrino flux times cross section integrated over At and 
AE, and is expressed as number of neutrino interactions per target: 

'S'Scut = {rilmin - /bkAt)/(Mact ' Np ■ e) 

The values of S are shown in table [H for the two LVD thresholds of the im- 
itation frequency, i.e., = 1 per 100 years and = 1 per month, two 
different masses, Mact = 1000 t and 330 t, and two values of -Ecut- As it can 
be seen an important improvement can be obtained by increasing the energy 
cut from 7 to 10 MeV. As was shown in figure [TJ the minimum cluster multi- 
plicity, for example at = 1 per 100 years, goes from 22 to 10, allowing an 
improvement of almost a factor of two in the sensitivity S-^^^^- 
It must be noted that in the on-line algorithm described so far we have ne- 
glected the capability of LVD to detect both products of the IBD reaction (see 
section 2.1). We can consider the signature of the reaction to build different 
burst selection algorithms. For example, we can require that all the H signals 
in the cluster are "signed", i.e., accompanied by a delayed L one (algorithm 
IBD-A). In the definition of the imitation frequency (eq. 1) fbk is substituted 
by fbk ■ PLbk (being P^b^ = 0.13 the probability for a H signal to be followed, 
in the same counter, by a L one due to background, averaged over the whole 
array) and, in (eq. 2), eE^^t by e^cut^n- The sensitivity of the algorithm IBD- 
A is shown in Table 1: even if the minimum multiplicity is lower, and the 
background rate is reduced, because of the n-capture efficiency (en = 0.5) the 
IBD-A method has comparable effectiveness or even less then the on-line one. 
We can also build several different algorithms, intermediate between the on- 
line and the IBD-A ones, requiring that only a fraction of the H signals in 
the cluster are accompanied by L ones (IBD-B)[I]. However, even if the IBD-B 
method efficiency results higher (as can be seen in Table 1, where an exam- 
ple is given), it does not exceed the on-line one enough to justify the loss of 
simplicity of the on-line algorithm and its independence from the model of 
supernova neutrino emission. Moreover, the on-line algorithm is sensitive to 
all possible neutrino interactions in LVD, both in the liquid scintillator and 
in the iron structure (that can represent up to 15% of the total number of 
interactions) . 



^ We reject all the clusters which have a number of "signed" H signals < k, such 
that: X]r=o -^(^' "^'P) — -^0; where P{r,m,p) is the binomial probability to have r 
signed pulses in a cluster of multiplicity m. We choose Pq = 0.1 as an example, the 
case Pq = corresponding to the on-line algorithm. 
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4 Discussion and conclusions 



The Large Volume Detector, located in INFN Gran Sasso National Laboratory 
has been designed to detect neutrino bursts from galactic gravitational col- 
lapses. In this paper we have described how the trigger and the event selection 
have been optimized to recognize on-line a neutrino burst from a supernova 
explosion even in the absence of an " external trigger" , such as the presence of 
an optical counterpart. The fast identification of such a neutrino signal is of 
utmost importance especially in view of the LVD participation to the SNEWS 
project, which should promptly alert the worldwide network of observatories 
to allow the study of the rare event since its onset. We have discussed the 
on-line algorithm currently in use at LVD and we have defined its sensitivity 
S in terms of minimum number of interactions per target unit requested to 
produce a burst alarm. 

If one assumes a model for the neutrino emission and propagation, it is pos- 
sible to express the sensitivity in terms of physical parameters of the source, 
such as its distance or the emitted neutrino flux. In particular we will use here 
the signal detected by Kamiokande-II and 1MB for SN1987A. Following [13] 
(see table 1 herein), [2] and [15], we adopt these values for the astrophysical 
parameters of SN1987A: average energy < E^^ >=14 MeV; total radiated 
energy E^j = 2.4-10^^ erg, assuming energy equipartition; distance D=52 kpc 
and average non-electron neutrino energy 10% higher than z/e[16j. Concerning 
neutrino oscillations (see [B] for a discussion), we consider normal mass hier- 
archy. We calculate the number of inverse beta decay signals expected from a 
SN1987A-like event occurring at different distances, for i^^ut = 7 and Ecut = 10 
MeV, and for two values of the detector active mass, Mact = 330 t and Mact = 
1000 t. For example, the number of detected IBD events for such a supernova 
in the center of the Galaxy (D = 10 kpc) is 230 for Ecut = 7 MeV and Mact = 
1000 t. Taking into account Poisson fluctuations in the signal multiplicity, we 
derive the on-line trigger efficiency as a function of the distance: this is shown 
in figure [9] (lower scale) for LVD working stand-alone and in the SNEWS. 
On the other hand, if we fix the distance of this supernova, e.g. at 10 kpc, we 
can derive the LVD sensitivity in terms of the minimum neutrino intensity at 
the source, varying the total emitted energy. The on-line trigger efficiency, as 
a function of neutrino luminosity (in terms of percentage of SN1987A one) is 
shown again in figure [9l upper scale. For example LVD working stand-alone, 
with 1000 if: of active mass and for Ecut = 10 MeV, is sensitive (at 90% c.l.) to a 
neutrino luminosity equivalent to 6% of that of SN1987A placed at a distance 
of 10 kpc. 
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Fig. 9. On-line trigger efficiency versus distance (lower scale) and per- 
centage of SN1987A signal at 10 kpc (upper scale) for £^11=7-10 MeV 
(light green and dark blue lines, respectively) and Mact = 330 t (dotted) 
and 1000 t (continuous). LVD stand-alone on the left panel and in the 
SNEWS on the right one. 



We can conclude that, without introducing any further check on the time 
structure, energy spectra and v flavor content of the signals in the cluster 
(which are postponed to the off-line analysis), LVD is able to identify on- 
line neutrino bursts from gravitational stellar collapses occurring in the whole 
Galaxy (D < 20 kpc) with efficiency > 90%. Such a sensitivity is preserved 
even if the detector is running with only one third of its total mass and stand- 
alone, with a severe noise rejection factor (<1 fake event every 100 years). The 
LVD on-line trigger efficiency can be extended up to 50 kpc (corresponding to 
the Large Magellanic Cloud) by introducing a cut on the visible energy at 10 
MeV. 
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Table 1 

Burst sensitivity for different selection algorithms, for two energy thresholds Ecut 
and for two values of imitation frequency; Mad is expressed in ton; rrimin represents 
the minimum cluster multiplicity (in parenthesis the minimum number of requested 
"signed" H signals for the IBD-B algorithm). 
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